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The  stress  concentration  of  an  oblate  spheroidal  inclusion  parallel 
to  the  stress  and  deformation  axis  is  obtained  by  using  Eshelby  theory. 

This  is  conplementary  to  our  previous  study,  in  vAiich  stress  COTvsentration 
is  ancdyzed  for  an  oblate  spheroidal  inclusion  nomal  to  the  stress  and 
deformation  axis.  Effects  of  the  elastic  stiffness  and  the  aspect  ratio 
of  the  inclusion  on  the  stress  concentration  are  examined  in  detail.  The 
internal  stresses  inside  the  inclusion  and  at  the  natrix- inclusion  boundary 
are  calculated  cojisidering  inhaiDgeraity  and  plaistic  defomation  effects. 
Stress  concentrations  at  graphite  flakes  and  nodules  in  cast  iron  are  dis- 


cussed. 

\ 
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STRESS  CONCENTRATION  DUE  TO  iW  OBIATE  SPHEROIIAL  INCLUSION 


1.  INTRODUCnON 

A good  deal  of  attention  has  been  given  to  stress  concentrations 
produced  by  holes  and  cavities  in  naterials.  However,  only  limited  efforts 
I have  been  devoted  to  effects  of  inlcusions,  inside  which  elastic  stiffnesses 

differ  from  those  of  the  surrounding  natrix.  In  our  previous  paper  [1]  (to 
be  referred  to  as  Paper  I),  stress  concentrations  in  naterials  bearing  an 
oblate  spheroidal  inclusion  were  studied  by  using  Eshelby  theory  [2].  Three 
cases  were  considered:  (1)  misfit  effect  due  to  an  isotropic  dilatatioral 
strain;  (2)  inhomogeneity  effect  due  to  a difference  in  elcistic  stiffnesses 
of  the  inclusion  and  matrix;  (3)  plastic  defomation  effect  due  to  the 
presence  of  a rvxi-defomable  inclvision  in  the  natrix,  viiich  is  plastically 
[ deformed.  The  influervces  of  the  «ispect  ratio  and  elastic  stiffness  of  the 


inclusion  were  evaluated.  In  aralyzing  the  inhomogeneity  and  plastic  defoma- 
tion  effects,  we  restricted  our  attention  to  the  case  of  the  oblate  spheroidal 
inclusion  having  its  broad  faces  nonial  to  tha  direction  of  the  exterral 
stress  and  plastic  deformation  (N  configuration).  In  thus  paper,  we  wish  to 
extend  our  calculation  to  thie  oblate  spharoidal  inclusion,  whuch  is  parallel 
to  tha  exterral  stress  and  plastic  defonration  axis  (P  configuration). 

Results  of  thie  calcualtion  for  tha  misfit  effect  were  unaffected  by  thus 
change  in  orientation,  but  those  for  two  othar  effects  were  quite  different. 
The  internal  stresses  inside  the  inclusion  and  at  tha  natrix- inclusion 
boundary  were  obt£uned  as  a fuiKrtion  of  tha  aspect  ratio  and  differences  in 
tha  elastic  moduli  of  tha  natrix  and  inclusion. 

I 

2.  INCLUSION 

Suppose  Xj  and  X^  be  a Cartesian  coordinate  system.  At  its  origin 

! 


is  an  oblate  spheroidal  inclusion  (cf.  Fig.  1),  of  v^ch  boundary  is  given 


as 


(1) 


The  aspect  ratio  of  the  inclusion  is  defined  by  k = a/c. 

The  matrix  and  inclusion  are  assumed  to  be  elastically  isotropic.  E 
and  E*  are  the  Young's  noduli,  y and  y*  are  the  shear  moduli,  and  v and  v* 
are  Poisson's  ratios  of  the  iratrix  and  inclusion,  respectively.  In  the 
following  ccilculatians , the  shear  modulus  ratio,  m,  defined  by  m = y*/y  , 
is  enployed  in  order  to  describe  effects  of  a differing  elastic  constant 
of  the  inclusion.  To  siiplify  the  presentation  of  results,  Poisson's 
ratios  of  matrix  and  inclusion  are  assumed  to  be  equal  to  1/3,  far  v^ch  m 
is  also  equal  to  E*/E.  This  is  recisonable  for  most  of  the  practical 
^plications . 


3.  INTERNAL  STRESSES  IN  AND  AROUND  THE  INCLUSI(»J 

Eshelby  [2]  developed  a method  for  obtaining  the  internal  stress  of  an 
ellipsoidal  inclusion  in  an  infinite  isotropic  matrix.  Ei^n  strain  inside 
the  inclusion  is  uniform.  Ife  also  gave  an  equation  to  evaluate  the  inhorao- 
geneity  effect  under  a uniform  external  stress  field  by  using  the  concept 
of  an  equivalent  inclusic»n.  As  in  Paper  I,  his  theory  is  used  her>e  to 
obtddn  the  internal  stresses  in  and  ^u^ound  the  inclusion. 

In  Paper  I,  only  tvro  independent  principal  stress  conponents  needed  to 
be  considered.  In  this  study,  however,  all  three  are  required  since  one  of 
the  geometrical  synmetries  is  lost.  This  ch^mge  adds  one  more  unknown 
oorponent  of  the  eigen  strain  e^^  of  the  equivalent  inclusion  in  comparison 


p 


to  the  previous  study.  After  solving  the  proper  simultaneous  equations  with 
T 

respect  to  e^^,  -Ow  principal  components  of  internal  stress  within  the  in- 
clusion, o^,  0^2  and  can  be  obtained  as  a function  of  the  aspect  ratio 
k and  the  shear  modulus  ratio  m = M*/y.  Intemca  stress  just  outside  of  the 
ixvclusion,  , can  also  be  obtained,  but  is  a function  of  the  position  at 
the  boundary.  We  will  present  primarily  the  results  of  at  the  equator 
point  Ag  along  axis  or  (0,  0,  c)  (cf.  Fig.  1). 

Two  effects  are  considered  in  this  study,  since  internal  stresses  due 
to  misfit  effect  can  be  obtcdned  for  the  present  geometry  ty  interclanging 
the  indices  1 and  3 in  the  results  presented  in  Paper  I: 

1)  Inhomogeneity  effect:  When  uniaxial  external  stress  is  applied 

A 

alcxvg  Xg  axis  at  infinity,  it  produces  elastic  strain,  ev^,  given  by 


(2) 


A A 

vdiere  e is  equal  to  a /E. 

2)  Plastic  defannation  effect:  When  the  inclusion  is  non-defomable 
plastically,  and  the  surrounding  matrix  undergoes  plastic  deforrration  along 


X-  axis  by  the  amount  of  e_»  intemeil  stresses  in  and  around  the  inclusion 
j P 

are  to  be  evaluated  as  a function  of  following  the  method  of  Tanaka  and 
Mori  [3],  In  this  c£ise,  the  eigen  strsdn  et.  inside  the  inclusion  is 


4.  RESULTS 


4.1  Inharpgeneity  Effect 

Results  of  intemcLL  stress  evaluation  due  to  inhonogeneity  effect  for  seven 
different  m values  are  shown  in  Figs.  2 to  4.  All  the  stresses  are  nor- 
malized by  the  applied  stress  d^.  When  inhoraogeneity  effect  is  cibsent 
(m  = 1),  vge  have  = 0 and  “ o regardless  of 

the  aspect  ratio  k,  as  expected.  As  k approaches  zero,  ^nd  (022^inh 

vanish  regardless  of  m,  vAiereas  approaches  an  asynptotic  values  of 

Figure  2 indicates  that  the  magnitude  of  stress  noraal  to  the  toroad  face 
of  the  inclusion  is  generally  small,  aixi  that  the  sign  of  changes  at 

k = 0.37.  For  elastically  hard  inclusion  with  m > 1,  beccmes  tensile 

far  thin  discs  under  tensile  loading  alcxig  axis.  For  inclusions,  with  m > 1, 
coitpressive  (tensile)  stress  develops  along  X2  axis  vdien  tensile  (coin>ressive) 
external  stress  is  e^lied  eilong  axis  (cf.  Fig.  3).  The  nagnitude  of  the 

corpressive  stress  depends  on  both  k and  m.  For  k = 0.07  and  m = 10,  the 

I A I 

ratio  reaches  -0.473.  Results  for  (C33^^j^j^  shown  in  Fig.  4 indi- 

cate no  unexpected  behavior.  The  approach  to  the  asynptotic  value  is,  however, 

nuch  slower  for  larger  m Veilues  than  for  snaller  ones. 

M 

For  inhomogeneity  effect  at  the  bour»daiy,  at  the  equator  point 

Ag  is  presented  in  Fig.  5.  The  stress  (as  norroilized  by  a ) is  positive  when 
m is  greater  than  unity  and  increases  with  decreeising  k or  increasing  m.  For  m 
less  unity,  and  o have  oppDsite  signs.  It  is  ijii)ortant  to  note  that 

an  elastically  hard  plate-li3ce  inclusion  parallel  to  the  tensile  stress  axis 
produce  significant  tensile  stress  at  the  edge  normal  to  the  plane  of  the 
inclusion.  On  the  other  hand,  the  tensile  stress  developed  at  the  edge  of 
an  elastically  soft  inclusion  under  oorpressive  loading  is  nuch  smaller  in 
magnitude  in  oonparison  to  a . 


4 


4.2  Plastic  Defornation  Effect 

Results  of  internal  stress  in  the  inclusion  due  to  plastic  deformation 
effect  are  shown  in  Figs.  6 to  8.  Here,  the  stresses  are  normalized 

by  E*ep,  v^iere  E*  is  the  Young's  modulus  of  the  inclusion  and  £p  is  the 
amount  of  plastic  defornation  along  the  axis.  For  a given  k,  a smaller 
m value  produces  larger  magnitude  of  cdl  principal  stresses  in  general. 

Both  (o?, ) and  (Ojo  ) have  the  opposite  sign  coirpared  to  e ; i.e.,  these 

Up  22  p " cr  IT  p 

stresses  are  conpressive  vdien  tensile  plastic  defornation  is  produced  in  the 
natrix.  Their  nagnitudes  are  snaller  for  larger  m and  approach  asymptotic 
values  of  zero  and  -0.19  with  vanishing  k,  respectively.  is  shown  in 

Fig.  8.  Its  nagnitude  is  greater  for  snaller  m and  approaches  the  asynptotic 


value  of  0.9375  E*ep  with  decreasing  k. 


For  plastic  defornation  effect  at  the  bouixlaiy,  "the  equator 

point  Ag  is  presented  in  Fig.  9,  Here,  the  stress  is  normalized  by  Ee^ 
where  E is  the  Young's  modulus  of  the  natrix.  The  effect  of  the  modulus 
ratio  m becomes  larger  at  snaller  k.  As  k goes  to  zero,  the  stress  approaches 
an  asynptotic  vadue  depending  on  m. 

4.3  0^  at  Equator  Points 

The  inclusion  considered  in  Paper  I has  the  synmetry  about  Xg  axis 
(also  the  stress  and  defornation  axis),  so  at  any  point  on  the  equator 
r^mins  unchanged.  However,  the  inclusion  considered  in  this  stutfy  does 
not  have  such  a synmetry.  So  far,  the  results  of  at  the  equator  point 
Aj  or  (0,  0,  c)  £u:«  presented.  It  is  necessary  to  exandne  the  variation  of 

0^^  CTi  the  equator  defined  by  = c^.  Consider  a unit  vector,  n, 

— ..22 
normal  to  the  equator,  where  n = (0,  n2,  n^),  satisfying  nj  + = 1. 

M 

According  to  Eshelby  [2],  on  the  equator  can  be  obtained  as 
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(4) 


M IT  1/T  2.T  2., 

■ 4 ^ ®kk^  ■ 2 ®)<k  ■ 2 ^ ®22  "2  ®33  *^3^^ 


M 2 

Dquation  (4)  can  be  written  as  ^3’ 


and 


S " 4 ^ * ®kk^ 


C = 1 E (e*^  - e"^  ) 

4 8 ^®22  33^* 


M 


1 , T ^ . Tv, 

2 ^*^22  *^33  ^ 


(4a) 


(4b) 


Since  0 £ £ 1,  has  a maxiram  at  either  n^  = 0 or  n^  = 1.  Here,  = 1 

corresponds  to  the  equator  point  A^,  and  n^  = 0 to  the  equator  point  A2  or 
(0,  c,  0)  (cf.  Fig.  1).  at  the  equator  point  is  then  equal  to  C^.  Values 
at  for  inhomogeneitv  and  plastic  defomation  effects  were  calculated. 
These  are  presented  against  k for  various  m values  in  Figs.  10  and  11. 

Figure  10  shows  the  results  for  inhomogeneity  effect.  The  stress  ratio 

is  always  negative  far  m > 1,  indicating  that  reaches 

a maxinum  at  the  equator  point  A^  under  tensile  applied  stress  (.a^  > 0). 

Under  coirpressional  ^plied  stress  (a^  < 0),  is  naxiraum  at  the 

equator  point  A^.  However,  its  nagnitude  is  nuch  snaller  than  that  at  A^ 
for  />0.  When  m is  less  than  unity  and  > 0,  at  A^  reaches 

a naxinum,  but  <^i>  at  ary  point  on  the  equator  is  alway  less  than  that 
of  . Therefore,  the  stress  concentratican  due  to  inhonDgerreity  effect  at 
the  edge  of  a soft  inclusion  parallel  to  the  stress  axis  is  deemed  insigni- 
ficant. 

Figure  11  shows  at  A2  arising  from  plastic  defomation  effect. 

When  these  results  are  oonpared  with  those  in  Fig.  9,  it  is  found  that  the 
location  of  a maximum  in  depends  on  the  sign  of  regardless  of  the 
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nodulus  ratio  m.  For  tensile  strain  (e  >0),  the  naxinum  occurs  at  A., 

P 

M 

whereas  reaches  a naxinum  at  for  compressive  strain  (e^  < 0). 

M 

Again,  the  nagnitude  of  at  A^  is  snaller  than  that  at  A^. 


4.4  CoBoparison  with  Previous  Results 

Results  obtained  here  for  a sphericcil  inclusion  agree  with  those 
reported  in  Paper  I.  For  this  special  case,  no  geometrical  distinction 
exists.  It  should  be  mentioned  here  that  Tana3ca  and  Mari  [3]  first 
evciluated  the  plastic  defornation  effect  quantitatively  using  Eshelty  theory. 
Originally,  they  considered  woric-Jardening  of  composite  materials.  Later, 
Tanaka,  Msri  and  NaJcanura  [4,5]  extended  the  calculations  to  deoohesion 
process  of  an  inclusion  in  naterial.  The  shape  of  the  inclusion  considered 
included  a thin  disc  parallel  to  the  deformation  axis  (with  k = 0) , a 
sphere,  and  a long  needle  along  the  X^  axis  (with  k = “). 

For  a spherical  inclusicm,  Tanaka  et  al.  [4]  obtained 

^'’aa^inh  “ ^ (7  - 5v)  + 03  - lOvJm  * 2(1  - 2v)  +'(1  + v)m^^  - v)  m c/^  (5) 
for  inhcitDgeneity  effect  and 


, I , (7  - 5v) .P*_ 

^‘^aa^p  ■ (7  - 5v)  (1  V)  + (8  - lOv)  (1  + vim  p 

for  plastic  deformation  effect,  respectively.  These  are  given  in  our 

I 

notation,  taking  v = v*.  Note  that  their  ejqjression  for  (o33)p  has  ca>e 
obvious  typographical  error,  vAiich  is  corrected  in  the  above  esquression. 

Our  results  are  in  complete  aigreeraent  with  theirs  given  in  Eqs.  (5)  and  (6). 
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For  an  infinitely  thin  disc  inclusion,  Tanaka  et  al.  [5]  obtained 


/_!  N - fi  ^ 2v  ^ A 


or  inhomogeneity  effect  and 


,J-\  _ (2  - V + V )E*e_ 


(7) 


(8) 


for  plastic  deformation  effect,  respectively,  again  taking  v = v*. 

Their  result  for  (033)p  agrees  with  ours  when  v is  taken  as  1/3.  However, 
the  result  for  ^033^^^^  diagrees  with  ours.  Apparently,  Eq.  (7)  is 
erroneous  because  the  equality  of  ^®33^£j^  = 0^  is  not  satisfied  even  vdien 
m = 1 unless  v = 0. 

5.  DISCUSSION 

It  is  known  that  graphite  flaJces  in  cast  iron  initiate  fracture  and  are 
responsible  for  the  low  tensile  strength  of  the  material.  In  most  previous 
studies,  the  graphite  was  treated  as  a cavity  (E*  = 0)  in  the  analysis  of  the 
stress  concentration.  However,  the  elastic  modulus  and  themHl  ejqsansion 
coefficient  of  graphite  are  obviously  not  zero.  Their  effects  on  the 
mechanical  strength  of  cast  iron  can  now  be  analyzed  more  ejqilicitly  by  using 
our  results. 

The  shape  of  graphite  flake  is  approximated  here  to  be  an  oblate  spheroid, 
vdiich  has  the  eispect  ratio  of  0.01  to  0.2.  The  limiting  case  of  nodular 
graphite  corresponds  to  k = 1.  TV#o  different  arrangements  for  the  disc  shaped 
gre^ihite  were  considered.  One  is  N configuration,  \diere  the  broad  face  of 
the  graphite  flake  is  nonnal  to  the  stress  axis  of  uniaxial  tension  or  ccxnpres- 
sion.  The  other  is  P configuration  where  the  graphite  flate  is  parallel  to  the 
stress  aucis.  These  are  referred  to  aa  ”N  graphite”  and  ”P  graphite”,  respectively. 
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The  maxiiium  stress  for  decohesion  of  graphite-iron  interface  is  expected  to 

be  stress  norrol  to  the  grpahite-iron  interface  at  the  polar  point  B or  rwrrral 

stress,  (parallel  to  the  inclusion  syTimetry  axis)  acting  at  the  equator  point  A. 

Three  different  effects  of  stress  concentrations  were  evaluated.  Young's 

moduli  of  graphite  and  iron  are  20.7  GPa  and  206  (3Pa,  respectively.  Thus, 

we  obtain  m = 0.1.  The  nagnitude  of  d^  was  taken  as  10”^E.  The  amount  of 

plastic  defomation  was  taken  as  1%.  It  is  assumed  that  the  natrix  ceases 

to  allow  plastic  relaxation  below  800°K.  The  thermal  expansion  coefficients 

of  graphite  and  iron  are  1.8  x 10”^/®K  and  12.0  x 10”^/°K,  respectively  [6]. 

T -3 

This  gives  the  themal  misift  strain  of  e = 5 x 10  . Results  of  Paper  I 

were  used  for  the  Ccilcualtion  of  misfit  effect.  For  each  configuration  of 

graphite,  the  toteil  stress  due  to  three  effects  was  calculated  as  a function 

of  k,  taking  tensile  and  ccrpressive  nature  of  a into  account.  The  results 

are  presented  in  Fig.  12,  vdiere  each  curve  is  labeled  for  the  configuration 

of  graphite  flake,  the  position  of  internal  stress  as  well  as  the  direction 

M I 

of  applied  stress.  Here,  (P,  A)  refers  to  at  A^  and  (P,  B)  to  at  B 

M I 

in  this  study,  vAiereas  (N,  A),  refers  to  0^3  ai  A and  (N,  B)  to  at  B in 
Paper  I.  The  subscript  T or  C refers  to  tensile  or  conprassive  applied 
stress.  Note  tlat  stresses  are  given  in  tenns  of  Young's  modulus  of  iron,  E. 

As  expected,  tensile  stress  at  the  edge  of  N graphite  (N,  A).p  is  quite 
high  urxler  tensile  loading  and  fracture  is  expected  to  initiate  at  the  edge 
regardless  of  k.  For  smaller  k,  ^ main  contribution,  but 

^'^3^inh  significant  as  k approaches  unity.  At  k = 0.01,  the  sum  of 

three  effects  is  20  x 10”^E,  which  is  five  times  smaller  than  ^ 

cavity  (m  = 0,  k = 0.01).  At  k = 0.1,  the  sum  is  40%  greater  than  ^0^3^;^^^ 
of  cetvity.  On  the  broad  face  of  N graphite,  (N,  B).p  is  still  tensile  for 
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A 

|h  k < 0.45  but  its  nagnitude  is  lower  than  that  of  a . Both  (N,  A)^  and  (N,  B)^ 

are  mostly  conpressive,  so  N graphite  is  expected  to  be  cin  inactive  source  for 
fracture  under  ccmpressive  loading. 

It  is  most  surprizing  that  both  (P,  A).j,  and  (P,  A)^  are  positive  and 
large  in  nagnitvide,  especially  for  k < 0.1.  This  iraolies  that  tensile  stress 
acts  at  the  edge  of  P eraohite  under  tensile  or  conpressive  loading.  Low 
tensile  fracture  strength  of  gray  cast  iron  is  thus  expected  frtxn  misfit 
effect  as  well  as  from  inhtxnogeneity  effect,  not  merely  frcm  the  latter.  Some 
plastic  defomation  in  the  vincinity  of  graphite  flakes  is  cilso  expected 
because  of  very  hi^  stress  level.  This  contradicts  our  assunption  and  in- 
validates a part  of  our  calculations.  However,  general  features  of  the 
results  are  still  applicable. 

Deleterious  effects  of  graphite  flakes  are  present  even  under  conpressive 
loading  according  to  the  present  calculation.  Still,  the  degree  of  stress 
concentration  is  nuch  less  uixier  oonpression  than  under  tension.  This  is  the 
obvious  origin  of  the  strength  differential,  effect  in  gray  cast  ir^  [7]. 
Another  obvious  consequence  is  the  favorable  effect  of  spheroidization  of 
graphite  [8],  as  is  well  known  industrially.  However,  the  stress  concentra- 
tion at  the  side  of  a nodular  grahite  is  more  than  50%  higher  than  ttet  of  an 
equivalent  cavity. 
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Fig.  1. 

Fig.  2. 

Fig.  3. 

Fig.  4. 

Fig.  5. 

Fig.  6. 

Fig.  7. 

Fig.  8. 

Fig.  9. 

Fig.  10. 
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Fig.  12. 


The  orientation  of  an  oblate  spheroidal  inclusion.  >^plied 
stress,  is  exerted  along  the  Xg  axis. 

Interral  stress  inside  the  inclusion  due  to  inhorogeneity  effect, 
against  k. 

Internal  stress  inside  the  inclusion  due  to  inhonogeneity  effect, 
^°'22^inh^^  against  k. 

Internal  stress  inside  the  inclusion  due  to  inhomogeneity  effect, 
^°L^inh''‘^  against  k. 

Internal  stress  at  the  equator  point  on  the  natrix- inclusion 
boundary  due  to  inhomogeneity  effect,  against  k. 

Internal  stress  inside  the  iiiclusion  due  to  plastic  defornation 
effect,  against  k. 

Internal  stress  inside  the  inclusion  due  to  plastic  defornation 
effect,  /E*c^  against  k. 

22  p p 

Internal  stress  inside  the  inclusion  due  to  plastic  defornation 
effect,  (022)p/E*ep  against  k. 

Internal  stress  at  equator  point  A^  on  the  natrix- inclusion 

M 

boundary  due  to  plastic  defornation  effect,  against  k. 

Internal  stress  at  the  equator  point  A^  on  the  matrix-inclusion 
boundary  due  to  inhomogeneity  effect,  against  k. 

Internal  stress  at  the  equator  point  A^  on  the  natrix- inclusion 
boundary  due  to  plastic  defamation  effect,  £p  against  k. 

Internal  stresses  at  equator  point  A and  the  polar  point  B of  iron 
natrix-grt^Jhite  inclusion  bourxlazy.  The  sum  of  the  thrae  effects 
(misfit  Cp  = 5 X 10"^,  inhomogeneity  m = 0.1  with  = lO”^  • E,  and 
plastic  deformation  Cp  * 1%)  is  shoMi  against  k for  two  oonfigurations 
of  graphite  inclusions  j N and  P greqihites.  Tensile  and  conpressive 
^lied  stresses  are  indicated  by  subscript  T aird  C,  respectively. 
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The  orictfitation  of  an  oblate  spheroidal  inclusion,  ^plied 

A 

stress,  o , is  exerted  along  the  axis. 


Fig.  2.  Internal  stress  inside  the  inclusic»i  due  to  inhonogeneity  effect 


efffec' 


Fig.  5.  Internal  stress  at  the  equator  point  on  the  natrix-inclusion 
boundary  due  to  inhonogeneity  effect,  against  k. 


effect,  (o:,,)  /E*e 


ide  the  inclusion  due  to  plastic  deformation 
against  k. 
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effect,  (Oot)  /Ee, 


due  to  inhOTDgeneity  effect,  (a,,)-  >,/cr  against  k. 


is  shown  against 
s;  N and  P graphites 


